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1. Introduction 

This paper completes a series of four reports on our calculation of Bk using HYP-smeared 
staggered fermions. In the previous reports we presented the results of fitting using SU(3) [|l|] and 
SU(2) staggered chiral perturbation theory, and our method for estimating the systematic error 
due to finite volume effects [^]. Here we focus on impact of the matching factors that we use 
to connect the lattice operators to their continuum counterparts, and explain how we estimate the 
errors that are introduced by truncating this matching at one-loop order. Further details will be 
given in Ref. [Q]. 



2. One-loop Matching 

To define matching factors we need to specify the continuum regularization and renormaliza- 
tion scheme used to define the operators in the continuum. If one matches perturbatively, as we do 
here, it is conventional to use MS regularization with the NDR (Naive Dimensional Regularization) 
prescription for 75 . One also needs to choose which class of lattice operators to use, and we follow 
earlier work and use the so-called two trace approach One then calculates the matrix elements 
of the Bk operator at some order (here one-loop) in both continuum and lattice regularizations. 
Equating them determines the matching factor — which is, in general, a matrix. 

Alternatively one can use the RI-MOM scheme, which is defined in any regularization, and 
determine the matching by a non-perturbative calculation on the lattice. The advantages and disad- 
vantages of this scheme are reviewed in Ref. ||^. We ultimately plan to use it, but so far have only 
used this method for bilinear operators 

Returning to the perturbative approach, the one-loop matching factors Z,y are defined through 

Of"-\^) = Ziji^,a)0Y'\a) (2.1) 

^ij = ^ij + ^ [7ij log (At«) + dj] (2-2) 
c,-,- =C^™'-CL.^" (2.3) 

where are continuum AS = 2 operators (here a single operator) renormalized at scale /i 

and OY^^ia) are the lattice operators required for the matching. Yij anomalous dimension 

matrix, while Cj^°"' and C|j-^" are the finite parts of the continuum and lattice matrix elements, 
respectively. The list of lattice operators which appear at one-loop order is given in Ref. [^]. This 
reference also calculates the matching factors for the HYP-smeared operators we use, but with the 
Wilson gauge action. The generalization to the Symanzik gauge action used to generate the MILC 
configurations will be presented in Ref.[§]. The results we use here are preliminary. 

In applying ( |2.1[ ) we make one simplification: from the rather long list of operators O^^'' which 
contribute at one-loop we keep only the four which have the same taste as the external kaons (1^5). 
This introduces 0{a/ (47r)) truncation errors which turn out to be of next-to-leading order (NLO) in 
SU(3) staggered chiral perturbation theory (SChPT) [[ig], and of NNLO in SU(2) SChPT 
Our SU(3) fits attempt to pick out the contribution of the missing operators and then remove them. 
Our SU(2) fits ignore these contributions as being of too high order. 
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Figure 1: Tree-level and one-loop matched Bk versus for degenerate valence quarks. Results are for the 
coarse (left) and fine (right) MILC lattices with m//m, = 0.01 /0.05. 



We now show how moving from tree-level to one-loop matching impacts the results for Bk- 
Here we first use what we call "parallel matching", in which the scale in the continuum operator 
is set to /X = l/a, and in which a, (which we take to be in the MS scheme) is also evaluated at 
this scale. The rationale for this choice is that it is a reasonable estimate for the typical momentum 
contributing in the matching. It is also possible to estimate the scale to use (usually called "q*") 
based on the integrand of the one-loop integrals, but we have not yet attempted this. 

In Fig. [l], we show the results for Bk on coarse (a ?a 0.12 fm) and fine (a w 0.09 fm) lattices 
before and after inclusion of the one-loop corrections. For clarity, we show only the points in 
which the valence quarks are degenerate. The roughly 20% reduction caused by the inclusion of 
one-loop contributions holds also for non-degenerate valence quarks. We also show the results of 
a four-parameter partial NNLO fit.^ 

The size of the one-loop shift (« 20%) is of the expected magnitude, given that a^(l/fl) is 
« 0.33 and 0.27 on the coarse and fine lattices, respectively. It should be kept in mind that, however, 
that Bk is scale dependent, and so the one-loop correction can be made larger or smaller by varying 
the scale chosen in the continuum operator. In other words, there is no precise way of defining the 
size of the correction. 

A noteworthy feature of these results is that the curvature at small Mk is larger after one-loop 
matching. This is even more pronounced in the results on the superfine lattices (a ^ 0.06 fm), shown 
in Fig. ^. Indeed, one can see that the fit to the tree-level results has an upwards "hook" at very 
small Mk, which is the result of the fit requiring a significant contribution from the taste-violating 
operators which are present because of truncation (and discretization) errors. These contributions 
behave as — logM/f in the chiral limit, and are finite there because the kaon that appears has non- 
Goldstone taste and so its mass does not vanish in the chiral limit. We expect such contributions 
to be of leading order in SChPT for tree-level matching, but of NLO, and thus much smaller, for 

'Details of this fit which will be explained in Ref. [B], and are not pertinent here. 
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Figure 2: Bk vs. M| for the MILC superfine lattice with m/ /m^ = 1/5. 



one-loop matching. This is consistent with our results. We stress that the curvature seen in the 
one-loop curves is not a surprise as the chiral logarithm has a fairly large coefficient. 

The subset of the data most relevant to extrapolating to the physical kaon is that in which 
the valence kaon is maximally non-degenerate. This effect of one-loop matching on this subset of 
the data is illustrated by Fig. ||, where we present the results of the SU(2) SChPT fits Bk on the 
coarse MILC lattices. We show an example of the "X-fit" (the extrapolation inM^ for fixed valence 
strange-quark mass) and the "Y-fit" (the extrapolation in the valence strange quark mass). These 
fits are explained in Ref. 



spacings, with taste-breaking discretization and truncation errors removed. In order to compare 
these values we next run them from 1/a to a common scale, which we take to be 2 GeV. Here, we 
use two-loop RG evolution 



3. RG Evolution 



The fitting procedures described in Refs. ||T], ^ result in values for (1/a) on the three lattice 




(3.1) 



Z 




where the anomalous dimension matrices y^'\ and the beta-function coefficients, Pi, are given, e.g.. 



in Ref. |11]. 



4 



Determination o/Bk using improved staggered fermions (IV) 



Jangho Kim 



o.e 



' ' ' ' I 

X tree 

O one loop 



0.7 



O.e 



0.5 



0.7 



0.6 



0.5 



' ' ' ' I 

X tree 

O one loop 



0.05 



0.1 



0.15 



0.2 



0.3 



0.35 



Xp (GeV^) 



0.4 

(GeV') 



0.45 



0.5 



Figure 3: Bk vs. Xp (left: X-fit) and vs. Yp (right: Y-fit) on the MILC coarse lattices (right: Y-fit) with 
mi/nis = 1/5. The [green] diamonds in the left plot show the result after extrapolation to the physical pion 
mass, Xp — M^, and removal of lattice artefacts from the fit function. In the right panel the [green] diamonds 
show the result of an extrapolation to the physical valence strange quark mass. 

After RG running, we have values of B/f (NDR,/x = 2 GeV) from the three lattices spacings, 
which still contain discretization and truncation errors. We attempt to remove the former by a 
linear extrapolation in a^, as described in Ref. [^]. As for the latter, we attempt to estimate these 
separately, as we now describe. 



4. Estimate of Two-Loop Terms 



Let B^^ be the value of Bk obtained using parallel matching (/x = l/a) at the fth loop level, 
and after extrapolation to the physical valence and sea-quark masses. Then we can define AB^^ as 



AB 



(i) 



B 



('•-1) 



B 



(i) 
K ■ 



(4.1) 



so that AB^^ represents the shift due to the /'th loop correction to Bk- We know fij^^ and and 
so we can calculate AB^^^; the results are collected in Table |I] for the SU(3) fits, and Table |2| for the 
SU(2) fits. One estimate of Afij^^ is then 

ABg)p.ABy'xa,(l/a), (4.2) 

with results also given in the Tables. 

We plot AB^' versus as{1-/a)'^ for the two analyses in Fig. ^ Linear fits yield intercepts 
consistent with zero. This is not surprising given that AB^^^ is obtained from a one-loop matching 



formula, eq. (2.1), in which the correction is proportional to as{l/a), and is then multiplied by 
a, ( 1 Z'^) again to obtain AB^ . The vanishing of the intercept is is not, however, an automatic result 
because the one-loop matching is applied before fitting and extrapolating the data, and involves 
contributions from lattice operators having different dependence on the quark masses. This means 
that AB^ need not be exactly linear in as{l/a). 
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Figure 4: AB^^^ vs. af for the SU(3) analysis (left) and for the SU(2) analysis (right). 
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An alternative estimate is simply to use 



(4.3) 



I.e. the naive estimate of the two-loop contribution. This turns out to be somewhat larger than 

(2) 

AB^ , as shown in the Tables. It also varies more rapidly with the lattice spacing. 

Since we extrapolate to the continuum limit assuming a linear dependence on a^, the more 
slowly varying truncation error does not extrapolate to zero, although it will be somewhat reduced 
from the value at our smallest lattice spacing. To be conservative however, we take the value of the 

(2)' 

truncation error /or the superfine lattices, and we use the larger of the two estimates, i.e. AB)^ . 
This gives the estimates that are included in the error budgets presented in Refs. |Q, ^]. 

There are various ways in which one can firm up and reduce the truncation error. One is to 
work on a yet finer lattice, which might allow one to fit to a combination of and a? errors, and 
will, in any case, reduce the size of the error. Another is to use two-loop matching. And, finally, 
one can remove all truncation errors, and replace them with statistical and some new systematic 
errors, by using non-perturbative renormalization. We are pursuing all three approaches. 



a (fm) 








as{l/a) 


<> 




0.12 


0.6898(59) 


0.5704(58) 


0.1194(83) 


0.3285 


0.039 


0.074 


0.09 


0.6118(95) 


0.5256(92) 


0.0862(132) 


0.2729 


0.024 


0.046 


0.06 


0.5963(83) 


0.5158(80) 


0.0805(115) 


0.2337 


0.019 


0.033 



Table 1: Results for tree-level and one-loop parallel-matched Bk using the SU(3) SChPT analysis (fit N- 
BT7). For all lattice spacings we use the MILC lattices for which mi jms =1/5. Also given are the one-loop 
shift, the value of aj, and estimates of the two-loop shift. 
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a (fm) 










ABf 




0.12 


0.6879(53) 


0.5751(49) 


0.1128(72) 


0.3285 


0.037 


0.074 


0.09 


0.6383(130) 


0.5358(122) 


0.1025(178) 


0.2729 


0.028 


0.048 


0.06 


0.5829(128) 


0.4937(119) 


0.0892(175) 


0.2337 


0.021 


0.032 



Table 2: As for Table |l| except using the SU(2) SChPT analysis (4X3Y-NNLO fit). 
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